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Abstract 

We study the scalar Higgs sector of the next-to-minimal supersymmetric stan- 
dard model with an extra U{1), which has two Higgs doublets and a Higgs singlet, 
in the light leptophobic Z' scenario where the extra neutral gauge boson Z' does 
not couple to charged leptons. In this model, we find that the sum of the squared 
coupling coefficients of the three neutral scalar Higgs bosons to ZZ, normalized by 
the corresponding SM coupling coefficient is noticeably smaller than unity, due to 
the effect of the extra U{1), for a reasonable parameter space of the model, whereas 
it is unity in the next-to-minimal supersymmetric standard model. Thus, these 
two models may be distinguished if the coupling coefficients of neutral scalar Higgs 
bosons to ZZ are measured at the future International Linear Collider by producing 
them via the Higgs-strahlung, ZZ fusion, and WW fusion processes. 
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I. Introduction 



Despite all its success against precision tests so far, the standard model (SM) is widely 
considered as the low-energy effective approximation of a fundamental theory. Two of the 
major directions to the fundamental theory beyond the SM may be via supersymmetry or 
grand unification, as numerous theoretical attempts to formulate it have been accumulated 
in these directions. 

The search for supersymmetric extensions of the SM is motivated by the No-Go the- 
orem, the gauge hierarchy problem, and the possibility of incorporating gravity by local 
SUSY [1]. Supersymmetric models extended from the SM require at least two Higgs dou- 
blets, in order to give masses to the up-quark sector and the down-quark sector separately. 
While the minimal version has just two Higgs doublets [2], nonminimal versions may have 
one or more Higgs singlets besides the two doublets. The next-to-minimal supersymmet- 
ric standard model (NMSSM) is the simplest one of those nonminimal versions, which has 
one additional Higgs singlet [3]. The NMSSM has five independent neutral Higgs fields, 
which may be separated into three scalar and two pseudoscalar Higgs bosons if CP is 
conserved. 

Grand unification is required in order to unify interactions in a consistent way, to 
reduce free parameters of the SM, and to explain naturally the fermion mass problem 
[4,5]. In 50(10), for example, the chiral fermion fields of the SM are unified into three 
generations of an irreducible representation. A number of Lie groups have been investi- 
gated within the context of grand unification. The mainstream of the grand unifications 
are such that their gauge groups are broken down, possibly several times, at higher scales, 
to yield the SM gauge symmetry as a subgroup. At the electroweak scale, the true char- 
acter of the grand unification may be investigated in retrospect if some structure of the 
unification group persist to remain. In other words, for instance, if the grand unification, 
after its breaking, gives rise to SU{2) x f/(l) x f/(l)' at the electroweak scale, the extra 
U{iy might lead the the original unification group. 

In recent days, supersymmetric grand unifications have also been studied in the liter- 
ature for the new physics beyond the SM. These studies have become more interesting by 
the recent numerical observation that, if one assumes the existence of superpartners of the 
SM particles, the unification of three renormalization-group gauge coupling coefficients 
works particularly well. The NMSSM may also be incorporated with grand unification 
such that at the electroweak scale the gauge symmetry of the NMSSM may accommodate 
an additional U{1) to become SU(2) x U{1) x f/(l)'. The string-inspired model might 
be a natural motivation to introduce the extra U{1) symmetry to the NMSSM. 

The NMSSM with the extra U{1) (UNMSSM) may differentiate itself from the NMSSM 
by exhibiting distinct phenomena [6,7]. The existence of Z' is a typical signal that distin- 
guishes the two models. Since the mass of Z' is generally expected to be larger than about 
600 GeV, establishing experimental evidences for Z' would take some time. Meanwhile, 
GzzSi — 1)2,3), the coupling coefficients of the neutral scalar Higgs bosons to ZZ, 
may provide a clue to distinguish the two models. In the NMSSM, there is a sum rule for 
these coupling coefficients [8]. We find that the corresponding sum rule is altered in the 
UNMSSM significantly. It is the purpose of this article to study the phenomenology of 
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the NMSSM with this extra U{1), which we may call the UNMSSM afterwards, and to 
discuss how to distinguish between them. 

At the future e+e~ International Linear Collider (ILC) with ^/s = 800 GeV, we antic- 
ipate that the neutral scalar Higgs bosons would most dominantly be produced through 
the WW fusion process, but also through the Higgs-strahlung process and the ZZ fusion 
process. We would like to show that we may tell, once the mass of at least one of the 
three neutral scalar Higgs bosons and its coupling coefficient to a ZZ pair are determined, 
whether the neutral scalar Higgs boson belong to the UNMSSM or to the NMSSM . 

II. UNMSSM 

We assume that the extra U{1) emerges from Eq. In other words, a low-energy superstring- 
inspired Eq model is assumed to give rise to SU{2) x U{1) x U{1)' of the UNMSSM. As 
Green and Schwarz has shown some years ago, string theories in ten dimensions can be 
anomaly free at the string-theory level if the gauge group is either Eg X Eg or SO{32) [5]. 
Between the two gauge groups. Eg x E^ can be reduced to Eq as an effective grand unifi- 
cation at low energy by compactifying additional six dimensions, and it can accommodate 
chiral fermions. 

The Eq may be decomposed into 

Eg D 50(10) X U{1)^ D SU{5) x U{1)^ x U{1)^ , (1) 

where SU{5) is further broken down to the SM gauge group, SU{3)c x SU{2)l x U{1)y- 
At the electroweak scale, the desired extra U{1) symmetry may be given as an orthogonal 
linear combination of U{1)^ and U{1)^ as 

f/(l)' = cos9eU{1)^ + sm9EU{l)^ . (2) 

The electroweak gauge group we consider is thus G = SU{3)c x SU{2)l x U{1)y x U{1)' . 
Sometimes, particular values of the mixing angle 9e designate specific model names as 
follows: the x-model for 9e = 0, the ip-model for 9e = 7r/2, the //-model for 9e = 

tan~'^(— y5/3), and the z/-model for 9e = tan~^ v^TS. We take the ?7-model. 

The structure of the neutral Higgs sector of the UNMSSM is identical to that of 
the NMSSM: There are two Higgs doublets Hi = {H^.H^) and H2 = {H^,H^), and a 
neutral Higgs singlet S. There are therefore ten real degrees of freedom. After spontaneous 
symmetry breaking, the three neutral Higgs fields develop vacuum expectation values as 
(Hi) = vi, {H2) = V2, and (5*) = s, where we assume that CP is conserved. For simplicity, 
we take only the third generation of quarks into account. Then, the superpotential for the 
Yukawa interactions in the UNMSSM for quarks and the exotic quarks may be expressed 
as 

W ^ htQ^eH^t'^R + hQ'^eHib% + hSkLkn + XH^eH^S , (3) 

where e is an antisymmetric 2x2 matrix with 612 = 1, and ht, hb and hk are respectively 
the dimensionless Yukawa coupling coefficients of top, bottom, and exotic quarks, tfj and 
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6^ are the right-handed top and bottom quark superfields, respectively, Q is the left- 
handed SU{2) doublet quark superfield of the third generation, and the right and left 
handed singlet exotic quark superfields are denoted respectively as kji and ki. Notice 
that there is a cubic interaction among the three Higgs multiplets with a dimensionless 
coupling coefficient A. 

The low-energy particle content of an grand unification model can be contained 
in the fundamental 27 representation of -Eg, where the SM matter particles occupy 15 
components, the two Higgs doublets occupy 4 components, and the Higgs singlet occupies 
2 components. The exotic quarks are introduced in order to complete the fundamen- 
tal 27 representation of E^. Thus, the remaining 6 components of the fundamental 27 
representaion are occupied by the exotic quarks. 

Let us denote the effective hypercharges of Hi, H2, and S, respectively, as Qi, Q2, 
and Q3. They are defined as 

{i = 1, 2, 3), where are the U{1)y hypercharges of the corresponding Higgs fields, 6 is 
the kinetic mixing parameter, and are given by 

Q"! = cos 9EQf + sin OeQ! , 

where Q^, Qf, and Qf are the hypercharges of U{1)', U{1)^, and U{1)^, respectively. 
These hypercharges satisfy Qi + O2 + = by virtue of the f/(l)' gauge invariance of 
the superpotential under U{iy. 

The kinetic mixing parameter 6 is given as 6 = gu/g'i, where g[ is the gauge 
coupling constant and Qu is a new gauge coupling coefficient arising from the mixing 
between gi and g[. Without mixing between gi and g[, gn as well as S would be zero and 
Qi would reduce to Q'^. At the electroweak scale, the analysis of renormalization group 
equation yields that g[{mz) ~ 0.46 [9]. It is observed that, since the coupling of Z' to 
a charged lepton pair is proportional to the effective hypercharge of the charged lepton, 
Z' would not couple to charged leptons if the effective hypercharges of charged leptons 
could be rendered zero [10]. This leptophobic condition for Z' is satisfied if 5 = 1/3. In 
our analysis, we assume that Z' is leptophobic. 

The extra neutral gauge boson Z' is essentially the core characteristic of the UNMSSM. 
The mixing between Z and Z' is described in terms of a mixing angle 0, which is given 
by 

where niz = (gj + g^Wf^ is the Z mass, m|, = 2g'^v'^{Ql cos^ (3 + Ql sin^ (3) + 2g[^s^Q\ 
is the Z' mass, and = \f gj + glg'i^'^iQi cos^ P — Q2 sin^ f3). Note that is not a free 
parameter but it depends on g'l, Qi {i = 1, 2, 3), tan/3 and s. 

Both the Z' mass and the mixing angle in the UNMSSM are strongly constrained 
by experiments. According to the latest Review of Particle Physics, the lower bound on 
the Z' mass in the rj model is about 745 GeV from pp direct search [11] and 619 GeV by 
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electroweak fit [12]. It is remarked that, if tlie Z' is leptopfiobic, tliat tlie experimental 
sensitivities are mucli weaker, and searclies for a Z' via hadronic decays at CDF are unable 
to rule out a Z' with quark couplings identical to those of Z in any mass region [13]. For 
instance, as Fig. 3 in Ref. [14] suggests, the Z' mass may be as small as 100 GeV. From 
a phenomenological point of view, we are interested in a light Z' . A light Z' appears 
when the vacuum expectation value of the Higgs singlet, s, is about 700 GeV. A general 
comment on s is that larger s leads to heavier Z' in the present model. 

It is constrained by precision measurements that in the UNMSSM should be smaller 
than 2-3 xlO~'^. This constraint on may also be loosened if the Z' is leptophobic such 
that the mixing angle might be as large as 0.06 [14]. Anyway, a large s satisfies the 
experimental constraint on in the leptophobic case. 

III. Numerical analysis 

Let us take a reasonable parameter space for numerical analysis. We take 1 < tan /3 < 30, 
< A < 0.85, and 700 < s < 2000 GeV, where tan/3 = ^72/^^1, and v = ^vl + vl = 
175 GeV. The U{iy gauge coupling constant is taken as g[ = 0.46, motivated by the 
gauge coupling unification. The SUSY breaking mass, msusY, arising from the soft SUSY 
breaking terms, is assumed to be within the range from 100 GeV to 1000 GeV. 

The three neutral scalar Higgs bosons, denoted as Si = 1, 2, 3), are obtained as the 
eigenstates of the 3x3 mass matrix for three neutral real Higgs fields. The upper bound 
on ms^, the mass of the lightest neutral scalar Higgs boson Si, at the one- loop level in 
the UNMSSM, including radiative corrections due to the quarks and squarks of the third 
generation, is given by [7] 

< A\;^ sin^ 2/3 + m| cos^ 2(3 

+ 2g[^v\Qi cos' [3 + Q2 sin' + Am|^ , (5) 

where A is the dimensionless coupling coefficient for the cubic Higgs interaction, and Am|^ 
is the radiative correction, which comes from the one-loop corrections due to the quarks 
and squarks of the third generation. It is a function of s. Among the relevant parameters, 
the upper bound on ms^ depends significantly on the values of A and Qj- 

Numerical calculation shows that the upper bounds on m^^ becomes maximum at 
tan /3 ~ 1.5 and then decreases as tan /5 increases. This is mainly because the effect of the 
Higgs singlet and the D-teim contribution due to the extra U{1) at the tree-level Higgs 
potential interfere each other. This behavior is different from the MSSM, where the upper 
bound on rris^ increases as tan (3 increases. We find that < 162 GeV in the parameter 
ranges of our choice in the leptophobic scenario of the UNMSSM, at the one-loop level. 

Now, let us study GzzSi — 1)2,3), the couphng coefficients of the neutral scalar 
Higgs bosons to ZZ . In the SM, the corresponding quantity is given as 

GzzH = 92m z/ cos 9w , 

where 6w is the weak mixing angle, and H is the neutral scalar Higgs boson of the SM. 
We may normalize GzzSi by Gzzh such that GzzSi = GzzSi/Gzzh- Explicitly, the 
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normalized coupling coefficients may be written as 

Gzzs. (0) = cos [30r,Cl + sin (302^Cl + -^^Cj , (6) 

^(JZZH 

where Oij = 1,2,3) are the elements of an orthogonal transformation matrix that 
diagonalizes the 3x3 mass matrix for the three neutral scalar Higgs bosons, and Cj 
(i = 1, 2, 3) are dimensionless parameters defined by 

Ci = cos 6'^ + sin 6*^ cos — cos 6w sin 6w sin , 

6*2 = cos 6*^ + sin 9^ cos + ^^^^ cos 9w sin 9^ sin , (7) 

92 

Cs = 9'iQz sin 9w sin . 

Notice that in the above expressions the Z-Z' mixing angle is present in Q, which is 
responsible for the interference effect of the extra gauge boson. 

If = 0, we would have Ci = C2 = 1 and C3 = 0, and GzzSi{^) — '^'^^ POu + sin pO-zh 
which are exactly the normalized coupling coefficients of Si to ZZ in the NMSSM. The 
mixing angle is indeed zero if Z' decouples from Z, or if the extra U{1) is absent in the 
UNMSSM. In the latter case, the UNMSSM reduces trivially to the NMSSM. 

Now, we would like to remark that GzzSri^)) the three normalized coupling coefficients 
in the NMSSM, satisfy a simple sum rule: 

EGL5,(0) = 1. (8) 

i=l 

On the other hand, GzzSii'P)-, the normalized coupling coefficients in the UNMSSM, would 
not satisfy the sum rule, since would not vanish in general. Because of the non-zero 0, 
we would have 




E ^Iz5, (0) = cos^ + sin^ PGt + \ ^l. (9) 



and thus non-trivial mixing between Z and Z' . Therefore, by examining the value of 
G'|25.(0), one may in principle detect the effect of the existence of Z' , and thus 
distinguish the UNMSSM from the NMSSM. 

For numerical analysis, we need to know the values of Oij. They are generated ran- 
domly by varying the relevant parameters within the allowed ranges, with the constraint 
that the matrix O is an orthogonal matrix. Then, we evaluate Y3i=iG%zsX^)- ■'■^ -^^8- 
1, we plot the result for Yh=i G%zs i^) against tan/?, for four different values for s: 700 
GeV (solid curve), 1000 GeV (dashed curve), 1500 GeV (dotted curve), and 2000 GeV 
(dotted-dashed curve). It is quite clear that ELi GlzsM) ^ ^ UNMSSM, for the 

given ranges of parameters. The difference is more vivid for larger tan (3 and for smaller 
s. 
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Figure 1: I]f=i Crizs' (0) against tan/5 in the UNMSSM, for four different values for s: 
700 GeV (solid curve), 1000 GeV (dashed curve), 1500 GeV (dotted curve), and 2000 
GeV (dotted-dashed curve). 



Please note that, for s = 700 GeV, we obtain \^\ < 0.03-0.037 and 644 < mz' < 650 
GeV, depending on tan j3. Our choice of the parameter values yields the ZZ' mixing and 
the Z' mass within their experimental constraints in the leptophobic case [14]. 

One may quantify the difference between the UNMSSM and the NMSSM by introduc- 
ing 

In Fig. 2, we plot Aj against tan/3 for s = 700 GeV, where parameter values are ran- 
domly varied within the allowed ranges. The points marked by star, circle, and cross are 
respectively Ai, A2, and A3. For most of the parameter space, Aj are about 7 %, which 
is quite recognizable. The plot of Fig. 2, as well as the result of Fig. 1, suggests that 
the neutral scalar Higgs sector of the UNMSSM may clearly be distinguished from that 
of the NMSSM. 
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Figure 2: Aj against tan /9 for s = 700 GeV in the UNMSSM, where parameter values are 
randomly varied within the allowed ranges. The points marked by star, circle, and cross 
are respectively Ai, A2, and A3. 



Let us now study productions of these neutral scalar Higgs bosons in e~^e~ collisions 
at the ILC, where the center of mass energy of the e~^e~ system is sufficiently high in 
order to produce them. The important channels to produce them are the Higgs-strahlung 
process e~^e~ — > ZSi, the WW fusion process e~^e~ — > Ue^eSi, and the ZZ fusion process 
e~^e~ — > e~^e~Si. We denote the production cross sections via each of these processes in 
the UNMSSM as af (</>), a^{(j)), and cxf (0), respectively. 

These production cross sections are related to the cross sections for the production 
of the SM Higgs boson via the corresponding process. For the Higgs-strahlung process, 
(7/^(0) in the UNMSSM is related to the production cross section for the SM Higgs boson 

af (0) = iGzzsX<P)/GzzH?crg, = Gl.sM^im ■ 

Likewise, for the ZZ fusion process. 
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TABLE 1: The masses, the normahzed couphng coefficients to ZZ^ and the production 
cross sections of the three neutral scalar Higgs bosons via Higgs-strahlung, ZZ fusion and 
WW fusion process in the UNMSSM, for tan/3 = 10, mg = 500 GeV, A = 0.35, A = 500, 
and m.p = 200 GeV at the ILC with ^/s = 800 GeV. Those quantities with = are the 
NMSSM values, except for (0), which are also the UNMSSM values: (0) = (0). 



i 


1 


2 


3 


ms, GeV 


140 


201 


295 




0.9018 


0.0162 


0.0108 


Glzsm 


0.9699 


0.0185 


0.0114 


(7^(0) fb 


62.68 


0.620 


0.103 


af(0)fb 


18.02 


0.297 


0.156 


af(0)fb 


19.38 


0.336 


0.166 


af(0)fb 


6.362 


0.060 


0.010 


af(0)fb 


6.843 


0.068 


0.011 



Similar relationship may be established between the production cross sections in the 
NMSSM and those in the SM: 

^f(0)=GLs.(0)asM (10) 
^f(0)=G|^5.(0)cT|M 

Note that, however, for the WW fusion process, the effect of the extra neutral gauge 
boson is absent and we may set = 0. Thus, G^ry/sX't^) = Gy^^rs.{0). Further, it holds 
that GwwSii^) /GwwH = GzzSii^), where Gwwh is the coupling coefficient of the SM 
Higgs boson to WW. Consequently, we have, for the WW fusion process, 

^nO)=Ghs.WaZ, 

both in the UNMSSM and in the NMSSM. 

At the ILC energy, the WW fusion process would be the most dominant production 
channel in e~^e~ collisions for the lightest neutral scalar Higgs boson. Once the existence 
of a neutral scalar Higgs boson is established at the ILC, its nature will be examined to 
determine whether it is the SM Higgs boson or the lightest neutral scalar Higgs bosons 
in the NMSSM or in the UNMSSM. If it turns out to be the latter one, the experiments 
would determine its mass rris^ and the coupling coefficient G'^^zsii^)- 

These experimentally measured values would then constrain the free parameters of 
each model. In terms of G'^^Sii^)^ may obtain o"f(0) and (Tf^(O), which are the 
production cross sections of 5*1 in the NMSSM. On the other hand, in terms of G'^^zs 
one may also obtain erf (0) and (Tf^(0), which are the production cross sections of Si in 
the UNMSSM. Such a situation is summarized in Table 1. 

Notice that all of the three neutral Higgs bosons are not so heavy for the parameter 
values we take and thus they are kinematically within the reach of the ILC. It is evident 
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from Table 1 that the production cross sections in the UNMSSM are predicted to be 
reasonably smaller than those in the NMSSM, due to the interference effect by the extra 
neutral gauge boson. Also note that Y^^^iG^zs^i^) ~ whereas Gl^s, (0) < 1- 
These numbers may then be compared with the corresponding data at the ILC. Practically, 
it would be difficult to measure all these quantities during the first stages of operations 
at the ILC. Nevertheless, by comparing G|^5^(0) with G\zsi{'P) Table 1, finding a 
distinction between the UNMSSM and the NMSSM would be possible. In particular, the 
sum rule for the coupling coefficients may easily be examined. Any discrepancy of the 
sum rule from unity would distinguish the UNMSSM from the NMSSM, by confirming in 
an indirect way the existence of the extra neutral gauge boson. 

IV. Conclusions 

The extra neutral gauge boson in the UNMSSM is expected to manifest its existence in the 
Higgs sector of the model. There are three neutral scalar Higgs bosons in the UNMSSM. 
The mass of the lightest scalar Higgs boson is constrained from above by, among others, 
the effective charges of the Higgs doublets and the Higgs singlet. The upper bound 

on it is about 162 GeV, at the one-loop level, which suggests that it can be discovered 
at the ILC. For certain parameter values of the UNMSSM, such as those in Table 1, it is 
kinematically possible that the ILC may produce all of the three neutral Higgs bosons if 
they exist. 

The existence of the extra f/(l)' gauge group induces the mixing effect between Z and 
Z' in the production processes of these neutral scalar Higgs bosons. It reduces the sum 
of the squared coupling coefficients of the three neutral scalar Higgs bosons to ZZ, in 
particular the coupling coefficient of the lightest neutral scalar Higgs boson, such that the 
distinction between the UNMSSM and the NMSSM is practicable. We suggest that if a 
neutral scalar Higgs boson is produced via the Higgs-strahlung, ZZ fusion, or WW fusion 
processes are measured at the ILC, one may be able to examine whether it is identified 
as the lightest one in the UNMSSM or in the NMSSM. 
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